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SUMMAEff 

Eelations were analyblcally oLtaineii for ttLe prediction of radial 
distrlLutiona of yelocity and temperature for fully developed laminar 
flow of gases and of liquid metals in tuLes with fluid properties 
variahle along the radiiis . The relations are applicahle to Loth heat- 
ing anil cooling of the fluid. By use of the relations for velocity and 
temperature distributions, relations were obtained among Nusselt number, 
friction parameter, and ratio of wall to bulk temperature. 

The Nusselt number and friction parameter were found to be inde- 
pendent of Eeynolds nvoaber and Prandtl ntimber. The effects of ratio of 
wall to bulk temperature on Nusselt number and friction parameter could 
be eliminated by evaluating the fluid properties at specified tempera- 
tures in the fluid. 


inteouugtion 

Most of the analyses of laminar flow and heat transfer in tubes, 
for instance the uniform-wall-temperature solution given in reference 1, 
have been carried out under the assumption that the fluid properties do 
not vary with temperature. Somewhat more work on variable fluid prop- 
erties has been done for the laminar boundary layer of a flat plate. 

(See, for instance, references 2 and 3). In reference 4, velocity 
distributions for laminar flow in tubes with variable viscosity and some 
approximate relations for temperature distributions and heat -transfer 
coefficients are given. 

In the analytical investigation reported herein, which was con- 
ducted in the NACA Lewis laboratory, both velocity and temperature dis- 
tributions are obtained for fully developed laminar flow in tubes of 
gases and of liquid metals with variable fluid properties . These 
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relations are applicable to botli heating and cooling of the fluid. By 
use of these Telocity anfl temperature distributions, Nusselt numbers 
anfl friction parameters are obtained. The results should be applicable 
to flovr at points -where the ratio cf distance from entrance to diameter 
of the -tube is on the order of one— -t\fentie-fch of the Eeynolds number or 
greater . 
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SYMBOIS 

The following symbols are used in the report: 

specific heat of fluid at constant pressure, (Btu)/(lb)(°F) 
Inside diameter of tube, (ft) 

exponent, Talue of -which depends on variation of -viscosity of 
fluid -wl-bh temper atiire 

friction factor, 2Tq/(pu.-^) 

friction factor -wl-th density evaluated at t^. 

friction factor -with density evaluated at wall temperature 

acceleration due to gravity, (32.2 ft/sec^) 

heat-transfer coefficient, " ‘^ 3 ^^ (Btu)/(sec) (sq. ft) (°F) 

■thermal conductivl-fcy of flxiid, (Btu) (ft)/(sec) (sq. ft)(°F) 

■bheimal conduct i-vlty of fluid evaluated at -tj,, 
(B-bu)(ft)/(eeo)(sq. ft)(°F) 

thermal conductivity of fluid evaluated at tQ, 

(Btu) (ft)/(sec)(sa ft)(°F) 

ITusselt number, hl)/k 

Nusselt number -with conductivity e-valuated at -t^ 

Nusselt number with conductivity e-val-uated at t^ 

Nusselt number -with conductivity e-valuated at tg 
Prandtl nuniber, Cpgj./k 
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q. 

rate of heat transfer toward tube center per unit 
(Btu)/(sec) (sq. ft)) 

area. 


rate of heat transfer at wall toward tube center per unit area 
(Btu)/(sec) (sq. ft)) 

a' 

dimensionless heat -transfer parameter, Q./o^ 



r 

radius, distance from tube center, (ft) 



^0 

inside tube radius, (ft) 



r' 

dimensionless radius parameter, t/vq 



Ee 

EeynoMs number, pujjD/|ji 



EOj 

Eeynolds number with density and -viscosity evaluated 

at tj. 

ESq 

Eeynolds number -with density and viscosity e-valuated 

at -wall 

t 

absolu-be temperature, (9 r) 



H 

bulk or avereige static temperature of fluid at cross 
of tube, (9 e) 

section 


temperature in fluid, 3;(tQ - t^) + -t^, (*^) 



^0 

absolute -wall temperature, (‘^) 


■ 

t' 

dimensionless temperature -parameter, t/tg 




dimensionless bulk-temperature parameter, t^/tQ 



u 

velocity parallel to asis of tube, (ft/sec) 




bulk or average velocity at cross section of tube, 

(ft/sec) 

u' 

dimensionless -velocity parameter, [i.Qu/(TQrQ) 




dimensionless bulk velocity parameter, 



X 

number 



z 

distance from tube entrance, (ft) 



\ 

dimensionless parameter, koto/(qo^o) 
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5 dimensionleBS friction parameter, TQrQ/(n{jn) 

5q dimensionleBS friction parameter -wltli Tiscosity e-vnluated at tg 

8^ dimensionless friction parameter -with, viscosity evalunted at t„ 

p atsolute viscosity of fluid, (it-seo/sq. ft) 

atsolute -viscosity of fluid evaluated at t^, (it-sec/sq. ft) 

Pq atsolute viscosity of fluid at -wall, (ib-sec/sq. ft) 

p mass density, (it-sec^/ft^) 

tuli: or" average density at cross section of tute, (it-sec^/ft^) 
pQ mass density of fluid at wall, (it-sec^/sq. ft^) 

1 shear stress in fluid, (ib/sg ft) 

Tq shear stress in fluid at wall, (ih/sq. ft) 

AITALTSIS 

In the present analysis, equations for the velocity and temperature 
distrihutions in a laminar s-fcream of gases or of liq.uid metals throu^ 
a tube are derived for the fully developed houndary layer. Ii*can. these 
eq,uations, relations for the ITusselt rnmiher and friction factor are 
obtained . 


Flow of Gases 

Assumptions . - The analysis for gases applies when the following 
assumptions hold: 

1. The veloci-ty u- at any given distance from the wall is inde- 
pendent of distance along the tube. This assumption holds for fully 
de-ve loped flow when the rate of density change in the direction along 
the -tube is small. The effect of density change becomes less important 
as "the velocity is increeised. Laminar flow can exist at these velocities 
if -fche tube diameter is small. 

2. The temperature difference tQ - t at any distance from the 
wall is independent of distance along -fche tube. This condition holds 
for fully developed temperat-ure distributions -with un i form longitudinal 
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distrlljutioii of heat transfer in the tube -when the fluid properties are 
constant and holds for Tariahle properties -when the heat -transfer rate 
varies sli^tly along the tube. 

3. Frictional heating effects can he neglected. These effects are 
small when the Mach number for the flow is small compared with unity; 
the Mach number is generally small when laminar flow throng a tube 
exists . 


4. The viscosity, the theimial conductivity, and the density vaiy 
with temperature; the specific heat and the Prandtl number are constant. 
The variations with temperature of the specific heat and the Prandtl 
number of gases are of a lower order of magnitude than the variation of 
viscosity, thermal conductivity, and density. 

5. Heat transfer by conduction in the gas in the direction along 
the tube can be neglected. This heat transfer is generally small com- 
pared with that transferred by convection along the tube except for very 
small velocities . 

Telocity distributions . - For laminar flow the shear stress is pro- 
portional to the velocity gradient, or 


T 



( 1 ) 


The relation for the variation of shear stress with radius for unaccel- 
erated flow (assumption l) is obtained by eg.uating the shear forces to 
the pressure forces acting on a cylinder of fluid of arbitrary radius 
a.nfl differential length (reference 5) . This relation gives 


T 



From eq.uations (l) and (2) 


T 




( 2 ) 


(3) 


In order to convert this eq^uation to dimensionless form, the follow- 
ing quantities are introduced: 


u' 


■"cFo 


(4) 


r» 


r 

^0 


(5) 
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Then 


An* M^O du 
dr* “ Tq dr 

Eq.uatlon (3) heccanes, "when "written in dimensionless form, 


r* 


M- du* 
^l 0 dr* 


( 6 ) 


(7) 


The Tariation of the "Tlscoslty of a gas or liquid metal "vri-th 
temperature can he shown from "vlacosity data to he given hy an equation 
of "bhe form 



Suhsti-buting this value for V-/\^q in equation (7) and separating 
"variahles results in 

du’ = - dr’ 

4 - 


( 8 ) 


Integration of this eq.uation between a point in the fltiid and the wall 
resiilts in 



In order to solve this eq,uation, the relation he"bween t* and r* must 
he obtained. This relation is obtained in "the following section. 


Temperature distributions . - For obtaining temperatiire distri- 
butions, a heat balance is fl3?st •wri'fcten for an annulus of fluid -with 
inside radius r, o"ufcside radius r + dr, and length dz. The heat 
entering "the outer cylinder of the annulus is 2jt(r + dr)(q^ + dq.)dz and 
the heat lea"ving "the inner cylinder is 2it r q. dz. The net heat trans- 
ferred into the annulus "through the cylinders is, If differentials of 
higher order than the first are neglected, 2tc dz(rdq. + q,dr) or 

2it [S(rq.)/Srjdr dz. The heat transferred throu^ the ends of the annulus 

by conduction is neglected (assumption 5). The heat picked up by fluid 
passing throu^ "the ann"ulua is 2it r dr pgu c^ dt or 

2it r dr dz pg u Cp St/^z. The heat balance can then be "written as 
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5(rq) St 

“ =1- sT 

Hie heat 'balance for. the -whole cross section of the tube is 
2it rg go = ^0^ Pb S T^b 


( 10 ) 


or 


Ih’om assumption 2, 


or 


2<lo = ^0 Pb S Cp 


^(•fcp - -fc) ^ Q 
Sz 

^^0 St ^ 

Sz Sz Sz 


(11) 


( 12 ) 


Diyiding eguation (lO) by eguation (ll) and using condition (12) and 
constant Cp (assumption 4) give 


1 ^(r« g/gp) _ _p_ 

2 u-|j 


(13) 


The static pressure is uniform, across "the "tube. From the perfect gas 
law, 

i ^ (14) 

Pb ^ 

Substituting eguation (14) in eguation (13) and Integrating produce 
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Aaotlier expression for q.’ is obtained as follows: The law for 

beat conduction is 


q. = k — 

^ dr 

E(iuation (16) can be written in dimensionless form as 

k dt' 


I* = A. 


ko dr' 


where 


X = 


^O'^O 


< 10^0 


(16) 


(17) 


(18) 


With the assumptions of constant Prandtl number and constant specific 
beat (assumption 4), tbe law for tbe variation witb temperature of tbe 
conductivity must be tbe same as tbe law for variation witb temperature 
of tbe viscosity, or 

^d 


^ = fJLV = 

Substitution of eq.uation (19) into equation (17) gives 

. dt' 
q:' = Xt'^ ^ 

Equating equation (15) to equation (20) gives 

r. 

t' u-^j' dr' 

Tbis equation can be integrated to give 

‘1 


(19) 


(20) 



1 - t' = - 

X 


r' ^b' 


r'u' 


dr' 1 dr' 


( 21 ) 


r' 


Tn order to solve equation (2l), t-^' , nt)', and u' must be 

obtained in terms of t' and r' . 
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From equations (23) and (25) 



( 22 ) 


(23) 


(24) 


f25) 


(26) 


Wien eq.uation (26) and tlie value for u* from equation (9) are sutsti- 
tuted in eq.uation (2l), there results 
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Eq.uation (27) glTes the relation "between t' and r’ for Tarlous 
Tallies of X for laminar flow of gases. The eq.uation can "be solved by 
the method of iteration; that is, an assumed relation between t' and 
r* is substituted, into the ri^t-hand side of the equation and new 
values are calculated by solving the equation. The new values are then 
substituted into the equation and the process is repeated. After the 
relation between t' and r' for each value of \ has been obtained 
from eq.uation (27), the relation between u' and r' can be obtained 
from equation (9) . 


Flow of Liquid Metals 

Assumptions . - Assumptions 1, 2, 3, and 5, which were used in the 
analysis of laminar flow of gases, are also used in the analysis of 
laminar flow of liquid metals In tubes. Assumption 4 is altered, how- 
ever, inasmuch as the variations with temperature of the density and the 
thermal conductivity for most liquid metals are small compared with the 
variations of Tlscosity. In the present analysis, it is assumed that 
only the viscosity varies with temperature. 

Yelocity distributions. - In the derivation of equation (9) for 
velocity distributions in gas flow, variable density and thermal con- 
ductivity did not enter the analysis. Equation (9) therefore applies 
also to liquid metals . 

Temperature distributions . - The analysis of temperature dlstrlbu- 
tlons for liquid metals is the same as that for gases except that p/p^ 
in equation (13) and k/ko ^ equation (17) are both set equal to 1. 

The equation for temperature distributions is, with these alterations. 


2154 




This eq.Tiation can "be solTed in the same -weLj as equation (27). The 
dimensionless hulk: temperature t^' for liquid metals is calculated 
from 





t 'u'r'dr' 


u’r 'dr ' 


rather than from equation (23), which was used for gases. 


(29) 


Nusselt Numbers and Friction Parameters 

Nusselt numbers . - The Nusselt number with the conductivity evalu- 
ated at the wall temperature is defined as 


where 



- ^0 
^ = to - tb 


(30) 


(31) 


From equations (30) and (31) and the definition of X in equation (18), 
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"“0 = xu'-'ti,'7 


( 52 ) 


■ffliere t^' is o'btained from eauationa (23), (9), and (27) for gases or 
from equations (28) ^ (9)^ sn<i (29) for liquid metals. 

Similarly tlie Fusselt number with the conductivity evaluated at the 
bulk temperature is 


Nu-b 


(33) 


Ml - %') 

An arbitrary temperature in the fluid t^. is defined by 


or 




= t^' + Ml - -bij') 


(35) 


■where x is some number. !Phe TIuseelt number -with the conductivity 

2 


evaluated at t^ is 


rruj. = 


Ml - S')(W' 


(36) 


Friction ■pfl'^'WTneters . - The friction parameter with the ■vi.ecoslty 
evaluated at the wall temperature is defined as 


So = 


'O^O 


0 = 


f37) 


or 


1 

° " “s' 


(33) 


where Ut,' is calculated from equations (25), (9), and (27) or (28). 
The friction parameter -with the viscosity evaluated at any temperature 
■bj- is , ‘then 

1 






( 39 ) 
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The convent ional friction factor -witli density evaluated at the •wall 
fQ can "be shown, from, the definitions of fQ, bQ, and Esq, to he 

given hy 


f 


0 


48o 

'Beo 


(40) 


or 

48^, 
f = — — 


(41) 


Conventional Telocity- and Temperature -Dis trihut ion Eai^ameters 

The conventional velocity-distrlhution parameter u/u-^ can he 
written as 


u ^ u' 
% "h’ 


(42) 


where u' is calculated from eq^uation (9) and u-jj' is calculated from 
equations (25), (9), and (27) for gases, or equations (25), (9), and (28) 
for liquid metals. 


Similarly the temperatvire -distribution paimieter (to-t)/(tQ-t|j) can 
he written as 




(43) 


where t' is calculated from equation (2?) for gases or equation (28) 
for liquid metals, and t^’ is calculated from equation (23) or (29). 


Special Case of Uniform Fluid Eroperties 

Equations (9) and (27) or (28) can he integrated for the special 
case where the fluid properties are constant hy setting the value of t' 
on the right-hand sides of these equations equal to 1. The results of 
these Integrations are 


and 


u' = |(1 - r'2) 


(44) 
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Also, Ujj' is found "by integration of eq^uatlon (25) to "be 

, 1 
= 4 


and, from, equation (29), 


11 1 

1 


From equations (32) and (47), tbe Hiisselt number is 

48 


Hu = 


11 


From equations (38) and (46) 

5 = 4 

and, from eq.uations (40) and (49), the friction factor is 

16 


f = 


Be 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


EISDI/TS AKD DISCUSSION 
Laminar Flow of Gases 

Yeloclty and temperature distributions. - Telocity and temperature 
distributions for fvilly deweloped laminar flow of gases in tubes with 
wariable fluid prox>erties are sbown in figures 1 and 2, respectively. 

The profiles were calculated from equations (42) and (43), where the 
various quantities in the equations were calculated as Indicated after 
the equations. The exponent d was found from viscosity data to have 
an average value of 0.68 for air and most common gases for temperatures 
between 0° and 2000° F. 

It is seen from figure 1 that heat addition to the gas sharpened 
the peak of the velocity profile in the central portion of the tube, 
whereas heat extraction causes a flattening of the profile at the center 
of the tube. These results agree qualitatively with the sketch given in 
reference 6 (p.lB6). Inasmuch as the viscosity of a gas increases with tem- 
perature, the viscosity near the wall for heat addition is greater than 
that near the center of the tube so that the gas is slowed down near the 
wall and accelerated in the center of the tube for heat addition to the 
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gaa. The opposite effect is produced for heat extraction from the gas. 
The variation of density across the tube apparently has a smaller effect 
on the velocity distrihutions than does the variation of viscosity. 

The changes in sharpness of the temperature profile at the tube 
center with heat addition and extraction, as indicated in figure 2, are 
in the same direction as those for velocity profiles because the thermal 
conductivity varies in the same way with temperature as does the viscos- 
ity. The shapes of velocity and temperature profiles for inmlnfl-r flow, 
however, differ somewhat in contrast to turbulent flow where velocity 
and temperature profiles are similar (reference 7). Althou^ equa- 
tions (l) and (16) have the same form, the shear stress varies linearly 
■vrith radius whereas the heat transfer per unit area does not, and hence 
the shape of the velocity and temperature profiles may be expected to 
differ. In turbulent flow the profiles are so flat that it makes little 
difference what assumptions are used for the variations of shear stress 
and heat transfer with radius. 

Nusselt numbers and friction parameters . - Wusselt numbers fmfl 
friction parameters are plotted as functions of tQ/t-jj in figures 3 »nri 
4, respectively. The Nusselt numbers were calculated from equa- 
tions (32), (33), and (36) and the friction parameters were calculated 
from equations (38) and (39). It is seen from the figures that the 
parameters vary considerably with to/ti^ when the fluid properties are 
evaluated at either the wall temperature or the bulk temperature. The 
effect of to/t-b on the parameters can, however, be nearly eliminated 
by evaluating the conductivity in the ITusselt number at 
^-0.27 ~ “0*27(tg-t^) + bg and by evaluating the viscosity in the fric- 
tion parameter at "to. 50 = 0*58(tQ-t^) + t^. The fluid properties in 

the Nusselt number must be evaluated at a different temperature than that 
at which the fluid properties in the friction parameter are evaluated. 
I nasm u ch as the velocity and temperatiire profiles are dissimilar and the 
bulk velocity and the bulk temperature are wei^ted differently (see 
equations (22) and (24)). 


Laminar Flow of Liquid Metals 

yelocity and temperature distributions. - The velocity and tempera- 
ture distributions for laminar flow of liquid metals shown in figures 5 
and 6 were calculated in the same way as were those for gases . The 
exponent d was found from the limited viscosity data available for 
liquid metals to have an avereige value of about -1.6 for lead, bismuth, 
and sodium potassium alloy at temperatures between 300° and 1500° F. 

The changes in velocity -prof ile shapes shown in figure 5 are 
opposite to those for laminar flovr of gases, because the variation 
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■with teanperature of the 'vlacoalty of liquid metals is opposite to that 
of gases . It is seen in figure 6 that the effect of to/t-jj on the 
shape of the temperature profile is sli^t for liq.uid metals, because 
the thermal conducti'vlty is considered constant, so that any change in 
the temperature distribution must be caused by variation of viscosity. 

Nusselt numbers and friction parameters. - Uusselt numbers 
friction parameters are plotted as functions of to/t^ in figures 7 and 
8, respectively. 

The curves in figure 7 show that the friction parameter can be made 
substantially in dependent of tg/t^ by evaluating the viscosity at 
to . 54 = 0*54(to-t-j3) + t^|. This temperature is close to the temperature 
used for evaluating the viscosity in the friction parameters for gases. 

The curve in figure 8 shows the variation of Husselt number with 
to/t^ for liquid metals . The effect of to/t-^, on the Nusselt numbers 
cannot be eliminated as in the case of friction parameters, inasmuch as 
the the r m al conductivity of liquid metals is assumed not to vai*y with 
temperature. The variation of Nusselt number with to/t^ is caused by 
the variation of viscosity with tenperature and can be nearly eliminated 
by multiplying the Niisselt nvmbers in figure 8 by the anpirical correction 

for liquids usually given in the literature (See, for 

exanple, reference 6, p. 190.) The fact that the Nusselt number for 
liquid metals varies with tQ/t^ indicates that Nusselt numbers for gases 

are probably affected by variation of viscosity as well as variation of 
conductivity. Evaluating the conductivity at a certain tenperature in 
order to eliminate the effect of tQ/t^ seems to be only a convenient way 
of presenting the results. 

CcmpariBon of Laminar and Turbulent Flow of Gases 

A comparison of laminar and t\irb\xlent flow of gases is given in 
figures 9 and 10. 

In figure 9 are shown calculated laminar Nusselt numbers from the 
present investigation and calculated tiirbulent Nusselt numbers for a 
Prandtl number of 1 from reference 7. The laminar Nusselt number is 
independent of the Reynolds and Prandtl numbers, whereas the turbulent 
Nusselt number is dependent on them. In order to eliminate the effects 
of tQ/i>| 2 , the fluid properties including the density in the Nusselt and 

Reynolds numbers are evaluated at the average of the wall and bulk tenper- 
atures to. 5 for turbulent heat transfer and at a tenperature between 
the btQk temperature and the tenperature at the center of the tube t_o.27 
for 1 ami -nn-r heat transfer. The fluid properties for heat transfer in the 
transition region Should evidently be evaluated at temperatures between 
these two. 
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In figure 10, tlie laminar friction factors were ceilculated from 
equation (41) and tlie turtulent friction factors were taken from refer- 
ence 7 . The laminar friction factor f is independent of Prandtl 
number. (The laminar friction parameter S is independent of hoth, 
Reynolds and Bcandtl numbers.) For friction factors, the fluid prop- 
erties including the density for hoth laminar and turbulent flow are 
evaluated at temperatures close to the average of the wall and bulk 
temperatures . 


SUMMARY OF RESULTS 

The following results were obtained from an analytical investi- 
gation of laminar flow in tubes : 

1. Variation of fluid properties across the tubes for laminar flow 
of gases caused a sharpening of the peak of the velocity and temperature 
profiles at the center of the tube for heat addition to the gas And 

a flattening of these profiles for heat extraction. 

2. Variation of fluid properties across the tube for laminar flow 
of liquid metals caused a flattening of the velocity profile at the cen- 
ter of the tube for heat addition to the liquid metal and a shasrpening 
of the profile for heat extraction. Variation of fluid properties had 
but a slight effect on the temperature profile for liquid metals. 

3. The effects of ratio of wall to bulk temperature on Russelt 
numbers anfl friction parameters for laminar flow with variable fluid 
properties could be eliminated by evaluating the fluid properties in 
these quantities at various temperatures in the fluid. 

4. The velocity and temperature profiles, as well as Russelt number 
and friction parameter, were indepeiadent of Reynolds and Prandtl number 
for fully developed laminar flow with fluid properties that varied 
across the tube. 


Rational Advisory Coramittee for Aeronautics, 

Lewis Fli^t Propulsion Laboratory, 

Cleveland, Ohio, February 1, 1951. 
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Figure 6. - Temperature distributions for laminar 
flow of liquid metals in tubes. 
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Figure 8. - Variation of Hueaelt number with ratio of wall to bulk tomperatura for laminar flow of liquid 

metals with constant thomal oonductlvity. 
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Figure 9. - Variation of Nusselt number with Reynolds nximber for flow of gases with 
fluid properties Including density In Nusselt number and Reynolds number evaluated 
at temperatures of t_o,27 “ ~0«27 (tQ-t^j) + t^ for laminar heat trjmsfer and at 
temperatures of tg.g - 0.5 (tQ-t^,) + t|, for turbxilent heat transfer at Prandtl 
number of 1 , 
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Figure 10. - Variation of friction factor with Reynolds number for flow of gases with 
fluid properties including density in friction factor and Reynolds number evaluated 
at temperatures of to,s 4 - 0.54 (tg-tt,) + t^, for laminar flow and at temperatures 
of tg.s - 0.5 (tg-tb) + tjj for turbulent flow. 
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